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Abstract 

Examples are given which prove the ICARUS detector quahty through relevant 
physics measurements. We study the fj, decay energy spectrum from a sample of 
stopping /i events acquired during the test run of the ICARUS T600 detector. This 
detector allows the spatial reconstruction of the events with fine granularity, hence, 
the precise measurement of the range and dE/dx of the with high sampling rate. 
This information is used to compute the calibration factors needed for the full 
calorimetric reconstruction of the events. The Michel p parameter is then measured 
by comparison of the experimental and Monte Carlo simulated fi decay spectra, 
obtaining p = 0.72 ± 0.06 (stat.) it 0.08 (syst.). The energy resolution for electrons 
below ~ 50 MeV is finally extracted from the simulated sample, obtaining {E^^^ — 
^Mc)/^MC = 11%/ [MeV] e 2%. 
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1 Introduction 



The study of muon decay has played in the past a major role for the understand- 
ing of weak interactions, being the only accessible purely leptonic process. Muon decay 
was first described in a model-independent way by Michel [1] , using the most general, lo- 
cal, derivative-free, lepton-number conserving, four fermion interaction. For unpolarized 
muons, the decay probability is given by: 

f p, .) ^ 1.= (3(1 - .) . f p(4. - 3, + 3.^1^ + 1/ W + (1) 

where iV is a normalization factor: x = is the reduced energy (ranging from me/ E-^^ 
to 1); Eg and nie are respectively the total energy and mass of the electron produced in 
the decay; E'max = ?7t.^/2 is the end-point of the spectrum; f{x) is the term accounting 
for the first order radiative corrections assuming a local V-A interaction [2]; finally, p 
and T] are the so-called Michel parameters, defined in terms of bilinear combinations of 
the coupling constants of the general four fermion interaction, and hence depending on 
the type of interaction governing the decay process. For the Standard Model (SM) V-A 
interaction, the parameters take the values psM = 0.75 and rjsM — 0. 

The V-A assumption has already been confirmed in muon decay with high preci- 
sion [3,4] by the determination of the whole set of Michel parameters and complementary 
measurements. Figure 1 (left) shows the theoretical shape of the fi decay spectrum for 
various values of the parameters p and rj. As shown in the figure, and expected from 
inspection of Equation 1, the shape of the spectrum is more sensitive to p, since rj is 
weighted by l/x and hence determines the shape at low energy. The radiative corrections, 
shown in Figure 1 (right), determine the shape of the spectrum near the end-point and 
therefore the value of p is very sensitive to them. It has been shown that the overall effect 
of the radiative corrections on the value of p is of the order of 6% [2]. 

The Michel parameter p has been measured in the past by several groups (see Ta- 
ble 1). Peoples [5], Sherwood [6] and Fryberger [7] have measured p in the late 60's using 
the high energy part of the p decay spectrum assuming the V-A value 77 = 0. Derenzo [8] 
has measured the lower part of the energy spectrum, and combined his data with the 
previous results (essentially Peoples' measurement) into a common two-parameter fit, to 
obtain p with a precision of about 0.4%, which is the most accurate existing measurement 
with no assumptions for the value of r). These results have been obtained using dedi- 
cated experiments involving data samples of typically several hundreds of thousands of 
events. More recently, the data from electron-positron colliders have been used to measure 
the Michel parameters of the purely leptonic r decay near the Zq resonance [9]. These 
measurements are based on the analysis of samples that typically include several tens of 
thousands of events. 

We present here a further measurement obtained with the ICARUS detector, during 
its test phase (2001). ICARUS is a project, proposed in 1985 [10], for the installation of a 
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Fig. 1. Left: Muon decay energy spectra for various values of the Michel parameters including 
the V-A values. Right: Effect of the first order radiative corrections on the muon decay spectrum 
(p = 0.75, T] = 0). 

large liquid Argon (LAr) time projection chamber (TPC) in the Gran Sasso Laboratory, 
Italy, for the study of neutrino physics and matter stability [11]. The physics poten- 
tial of this type of detector has been extensively described elsewhere, both for the final 
project [10] and its initial phase with the ~ 600 t LAr prototype (ICARUS T600) [12]. 



The sample of events in which a muon enters, stops and eventually decays in the 
detector's sensitive volume -hereafter called stopping muon sample- constitutes an impor- 
tant benchmark to evaluate the physics performance of ICARUS. Because of their simple 
topology, stopping muon events are relatively easy to reconstruct in space, allowing the 
computation of the different calibration factors needed in the full calorimetric reconstruc- 
tion. Thus, we can study the muon decay spectrum and measure the Michel p parameter, 
which constitutes the first physics measurement performed with the novel ICARUS de- 
tector technology, and proves that the technique is mature enough to produce physics 
results. Our new result is not competitive with those obtained from n decay, and barely 
with those obtained from r decay. However, it must be remarked that this result has been 
obtained using 1858 muon decay events with a non optimized experiment. Our result 
stresses the capabihties of the ICARUS technology to produce robust physics results. 



2 Experimental setup 



ICARUS T600 [13] is a large cryostat divided in two identical, adjacent half-modules 
of internal dimensions 3.6 x 3.9 x 19.9 m^, each containing more than 300 t of LAr. Each 
half-module houses an internal detector composed of two TPC's (referred to as Left and 
iZi^/i^ chambers) , the field shaping system (race track electrodes), monitors, probes, photo- 
multipliers, and is externally surrounded by a set of thermal insulation layers. Each TPC is 
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Author 


Value 


Assumption 


Peoples 


0.750 




0.003 


r) = 


Sherwood 


0.760 




0.009 


r) = 


Prybcrgcr 


0.762 




0.008 


r] = 


Derenzo 


0.752 




0.003 


-0.13 <?] < 0.07 


SLD 


0.72 




0.09 ± 0.03 


lepton univers. 


CLEO 


0.747 




0.010 ± 0.006 


lepton univers. 


ARGUS 


0.731 




0.031 


lepton univers. 


L3 


0.72 




0.04 ± 0.02 


lepton univers. 


OPAL 


0.78 




0.03 ± 0.02 


lepton univers. 


DELPHI 


0.78 




0.02 ± 0.02 


lepton univers. 


ALEPH 


0.742 




0.016 


lepton univers. 


This analysis 


0.72 




0.06± 0.08 


-0.020 < r? < 0.006 



Table 1 

Results from previous measurements of the Michel p parameter. First (second) quoted error is 
of statistical (systematical) origin. Single error bounds correspond to statistical and systematic 
errors added in quadrature. 

formed by three parallel planes of wires, 3 mm apart, oriented at 0, ±60° to the horizontal, 
with 3 mm pitch wires, positioned along the left and right side walls of the half-module. 
The cathode plane is parallel and equidistant to the wire planes of each TPC A high 
voltage system produces a uniform electric field (500 V/cm) perpendicular to the wire 
planes, forcing the drift of the ionization electrons (the maximal drift path is 1.5 m). The 
electric field is uniform in the volume contained between the wire planes, the cathode and 
the planes of race track electrodes (hereafter referred to as LAr active volume) . 

The ionization electrons produced in the LAr active volume drift perpendicularly 
to the wire planes due to the applied electric field, inducing a signal (hit) on the wires 
near which they are drifting while approaching the different wire planes. By appropriate 
biasing, the first set of planes can be made non-destructive {Induction planes), so that the 
charge is finally collected in the last plane {Collection plane). Each wire plane provides 
a two-dimensional projection {view) of the event, where the position in one coordinate is 
constrained by the hit wire, while the signal timing with respect to the trigger gives the 
position along the drift direction. Each track is sampled by a large number of wires. This 
is one of the main feature of the ICARUS technology, exploited in the present analysis. 

A full test of the T600 experimental set-up on the surface of the earth was carried 
out in Pavia (Italy) during the period April- August 2001. One T600 half-module was fully 
instrumented to allow a complete test under real experimental conditions. Although the 
test was mainly intended as a technical run, a substantial amount of cosmic ray data was 
acquired, since the detector was not shielded against cosmic rays as will be the case in 
the Gran Sasso underground laboratory. The detector's conditions during data acquisition 
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Fig. 2. Run 966 Event 8 Right chamber: muon decay event views corresponding to the Collection 
(left) and second Induction (right) wire planes. 

(DAQ) were not stable, being subject to the optimization of various working parameters 
of the detector. 

One of the main goals of the test run was the detection of long muon tracks crossing 
the detector at large zenith angles (horizontal muons), extremely useful for an overall test 
of the detector performance. Therefore, an external trigger system was set up to select 
the horizontal muon events out of the overwhelming background of atmospheric showers. 
The system consisted of two plastic scintillator layers suitably positioned and arranged 
in a proper coincidence trigger logic. In spite of the optimization of the trigger for the 
acquisition of horizontal muon tracks, other kind of events (in particular, stopping muon 
events) are expected to be acquired within the DAQ time window. Most of these events 
can not be directly correlated to the triggering event, hereafter they will be referred to as 
out-of-time events (conversely, those which are correlated to the triggering event will be 
referred to as in-time events). 



3 Data selection 



Data selection was carried out by visual scanning using topological criteria. Stopping 
muons can be recognized as minimum ionizing tracks entering the detector, with increasing 
energy deposition and multiple scattering angle when approaching the stopping point. 
Two different topologies are expected depending on the process undergone by the stopped 
muon, namely: decay or nuclear absorption. The probability of the absorption processes 
has been evaluated to be 73% (0%) for negative (positive) muons [14]. For absorption 
events (/i~ -\- Ar CI* + u^), the excited CI nucleus decays emitting photons which then 
may interact via Compton scattering. For decay events, a minimum ionizing electron track 
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(length shorter^ than ~23 cm) follows the muon. In the present analysis, both decay 
and absorption events have been used in the determination of the calibration factors 
for the calorimetric reconstruction (see section 4.2). Figure 2 shows the two-dimensional 
projections (produced by the Collection and the second Induction planes respectively) of 
a typical muon decay event. 

In the case of decay events, the muon and electron tracks are separated by the muon 
stopping point, which is taken as the hit with the maximal ionization. The stopping point 
is assigned to the muon track. This may produce, for events where the muon and electron 
tracks partially overlap, a loss of the first hits from the electron track. This effect has been 
evaluated and introduced in the analysis using a sample of simulated events, as explained 
in section 5.1. 

A total of 5830 triggers were scanned, containing 4548 stopping muon events. Out of 
them, 3370 (74%) were further selected by the preliminary quality cuts. Among the 1178 
rejections, 45% of the cases were events with deficiently measured drift charge since the 
event happened in a region of the detector where the wires had an inadequate polarization 
during the run; 25% of the cases were due to an event occurring in a region with a 
substantial component of high amplitude, correlated noise, which may fake the presence 
of an ionizing track, hence, distorting the real event; the remaining 30% of rejections 
were due to several other effects, mainly events not fully contained in the LAr sensitive 
volume ^ , failures of the reconstruction program due to complicated topologies or absence 
of one or more views due to problems of the DAQ system. 



4 Data reconstruction 



4-1 Spatial reconstruction 



A full 3D reconstruction of the selected muon events is first performed. The spatial 
reconstruction of the muon tracks is needed in order to compute the calibration factors 
entering the calorimetric reconstruction of the events, as explained in section 4.2. 

A detailed description of the spatial reconstruction tools has been reported else- 
where [13,15]. The 3D reconstruction is performed in a hit basis using a three-step pro- 
cedure: 

(1) First, hits are searched for independently in every wire of each wire plane. Hits are 

^ The minimum observable track length is given by the wire pitch: 3 mm. 
^ An event is considered not to be fully contained if the muon or electron end points (identified 
by the raise in the energy deposition per unit length) are not present. No further fiducial cut on 
the LAr sensitive volume is required. 
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Fig. 3. Run 966 Event 8 Right chamber: fuhy reconstructed muon decay event. 



identified as signal regions of a certain width with output values above the local mean. 
No information from adjacent wires is used at this stage. A precise determination 
of the hit position and charge are carried out by means of a fit using an analytical 
function. 

(2) The identified hits are associated in the second step into 2-dimensional clusters of hits 
belonging to a common charge deposition, such as tracks or showers. Clusters provide 
criteria for the identification of the different patterns and for the discrimination 
between signal and noise hits, based on the cluster hit multiplicity. 

(3) In the third step, the 3D coordinates of the hits from the different reconstructed 
clusters are computed. Each wire plane constrains two spatial degrees of freedom of 
the hits, one common to all the wire planes (the drift coordinate) and one specific for 
each plane (the wire coordinate). The redundancy on the drift coordinate allows the 
association of hits from different views to a common energy deposition, and together 
with the wire coordinates from at least two planes, the determination of the hit 
spatial coordinates. 



Figure 3 shows a 3D reconstructed stopping muon event. The muon decay topology 
is visible, with three reconstructed tracks corresponding to the muon, the decay electron 
and a 7 — e conversion, respectively. 
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Fig. 4. Kinetic energy vs. range for typical well (top) and badly (bottom) reconstructed muon 
tracks. 



Calorimetric reconstruction 



The ionization charge is precisely measured at the Collection wire plane. The energy 
associated to a given hit is related to the collected charge by means of 

E=^ e(*-*")/- Q (2) 
R 

where C = (152 ±2) x 10"^ fC/(ADCx/xs) is the calibration factor [16]; W = 23.6l[]J eV 
is the average energy needed for the creation of an electron-ion pair [17]; R is the electron- 
ion recombination factor; (t — to) is the drift time of the electrons; Tg is the drift electron 
lifetime, which parametrizes the attachment of drift electrons to impurities in LAr; and 
Q is the measured charge. R, and are extracted from the reconstructed muon tracks, 
essentially by tuning them so that the measured energy corresponds to the theoretical 
expectation for stopping muons. This method determines the electron energy in a bias- 
free way, since all the calibration parameters are tuned using exclusively muon tracks. For 
the real experimental conditions to is given by the triggering system, but is unknown for 
most of the stopping muon events used in this analysis, since they are out-of-time events. 



Only those tracks successfully reconstructed in space are used in the determination 
of Te and R. A quantitative measurement of the quality of the spatial reconstruction is 
provided by the goodness of the fit of the muon track to the theoretical kinetic energy 
vs. range curve, when R, to and Tg are left as free parameters. A total of 2690 stopping 
muon events are selected using this method. Two examples of the energy vs. range curve 
measured for typical well and badly reconstructed muon tracks are shown in Figure 4. 
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The drift electron lifetime, Tg, is measured for data taking intervals of 24 hours. The 
measurement method is based on the observation of the charge attenuation with the drift 
distance, and has been published elsewhere [16]. However, two main differences must be 
stressed [15]. First, the charge of a hit depends on the particle momentum, but also on 
the angle between the track and the wire. Therefore, in general the hit charges are not 
directly comparable. In order to solve this problem, we normalize each charge with the 
theoretically expected energy, computed with the Bethe-Bloch formula given the range and 
length associated to the hit. Second, to is unknown for most of the events, and hence also 
the absolute drift time/distance. In order to merge hits from different tracks with different 
to into a common measurement of Tg, we use the relative charge attenuation (rather than 
absolute) as a function of the drift time/distance. Using this method, we obtain values of 
Tg ranging from 1.20 to 1.70 ms (depending on the data taking period), in agreement with 
the results previously published [16]. Data are also compatible with a slight dependence of 
the electron hfetime on the height: the electron hfetime is about 15% lower at the top than 
at the bottom of the liquid argon volume, with an almost linear dependence on height. 
The measured charge is corrected for the attenuation using the value of Tg averaged for 
the LAr volume, therefore this decrease is considered as an uncertainty when computing 
R and introduces a contribution to the total energy resolution ^ 

R is computed as a function of the theoretical dE/dx (estimated from the range using 
the Bethe-Bloch formula) , by comparing the measured charge (corrected for the finite drift 
electron lifetime) with the theoretical expectation for stopping muon tracks [15]. Using a 
sub-sample of 112 in-time events, for which to is known, we obtain -Rmip = 0.640±0.013 for 
minimum ionizing particles, where the error is dominated by the uncertainties of Tg and 
of the length of the track segments used to evaluate dE/dx. i^mip is used as a reference to 
tune the value of to for out-of-time events. The uncertainty of this procedure is estimated 
by applying it to in-time events: the values of to obtained for these events are Gaussian 
distributed, with mean 2 ±9 /is and width 82 ±8 fis. This can be interpreted as a negligible 
shift on the energy scale (< 0.3%) and an extra contribution to the energy resolution of 
^ ^ = 7%. However, it must be stressed that this contribution is due exclusively to 

the fact that the analyzed events are out-of-time, therefore it will not be present for the 
detector's real experimental conditions. The whole muon data sample can be used next 
to compute R for higher values of dE/dx up to ~ 5 MeV/cm. A linear dependence of R~^ 
on dE/dx with slope 0.11 ± 0.01 cm/MeV is found. 



5 Results 



5. 1 Determination of the Michel p parameter 



Prom the selected 2690 stopping muon events, a total of 1858 contain electrons 
available for the determination of the Michel p parameter (while the remaining ones are 
mainly absorption events) . The measured energy does not correspond directly to the p de- 
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Fig. 5. Energy spectrum of electrons from muon decays in the ICARUS T600 detector (electron 
energy lost by ionization in LAr). The plot shows the measured distribution (dots) compared 
both to the SM expectation (solid line) and the best fit with p and r] as free parameters (dashed 
line) . 

cay spectrum (Equation 1) since no attempt to recover the energy loss by Bremsstrahlung 
radiation has been carried out ^ . The spectrum corresponds to the fraction of the electron 
energy lost by ionization in LAr. Therefore, in order to measure the Michel p parameter 
we have chosen an approach based on the comparison with a Monte Carlo (MC) simulated 
event sample. This sample has been generated using FLUKA [18] and consists of 10 000 
electrons from muon decay events inside the detector's sensitive volume. The simulation 
includes all detector effects except for the presence of impurities in LAr. MC events are 
reconstructed using the same tools as for the data, and the effect of the impurities and the 
determination of are included in average by smearing the measured energy {E) using 
a Gaussian function of width a = 0.09E. The effect of the muon track (which is absent 
in the MC sample) is the loss of some of the first hits of the electron track due to the 
overlapping with the muon track. The number of hits of MC electron tracks is found to be 
on average 2.14 higher than for data electron tracks. Thus, the effect of the muon track 
can be included on average in the simulation by removing, when computing the energy, 

^ This does not decrease the statistical accuracy of the measurement and avoids the problem of 
associating photons to the electron track. Indeed, background conditions in data taking at the 
surface are such that it was preferred not to worry about a selection of photons pertaining to 
the electron. 
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Fig. 6. Results from the fit of the Michel electron spectra. The dot and band show, respectively, 
the position of the minimum and la region within the considered parameter's space. 

the first 2 hits (3 in 14% of the cases) of the electron track. The measured and simulated 
energy spectra are compared in Figure 5, and they are found to be in good agreement 
ixVndf = 14.0/20). 

Given the greater sensitivity of the shape of the spectrum to p than to rj, and the 
level of correlation between these two parameters (see Equation 1), we measure p while 
constraining rj within its experimentally allowed interval (—0.020 < 7] < 0.006) [19]. Since 
the MC sample is generated using the SM values (psM and 77sm), the spectrum for an 
arbitrary pair of values, p and r], is built by weighting each MC event with a factor: 

^ = dpf 7 (3) 

-^\Xmc] Psm, Vsm) 

where xmc = -SMcZ-E'max, and Eyic is the generated energy. We extract the value of 
p as that for which the best fit between the simulated and measured energy spectra is 
obtained, which yields p = 0.72 ±0.06, where the error has a statistical origin and includes 
the correlation with rj. Figure 6 shows the best fit point and the la region of the fit within 
the considered Michel p-rj plane. As expected, our measurement has no sensitivity for the 
determination of rj within the considered experimental bounds, and yields a stable value 
of p for the whole allowed r] range. 

There are two types of systematic uncertainties that can affect this measurement, 
namely: a bad estimate of the energy resolution and a systematic shift in the global energy 
scale. Using the MC sample, the uncertainty of p is estimated to be 0.01 (0.03) for an extra 
±5% (±10%) contribution to the energy resolution, and 0.04 (0.08) for a shift in the energy 
scale of ±1% (±2%). The energy resolution is determined essentially by the wire's signal- 
to-noise (S/N) ratio and the additional contributions from the drift electron lifetime, 
recombination and to corrections. All these contributions have been extracted from the 
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data and included in the MC sample, and its uncertainty is estimated to be smaller than 
5%. Furthermore, the effect on the final result becomes negligible by removing the last two 
bins of the electron spectrum when performing the fit. Conversely, the absolute energy 
scale is given by R, which is determined for minimum ionizing particles with an accuracy 
of 2% including systematic cfi^ccts (see section 4.2). Therefore, the total error is dominated 
at this level by the systematics due to the uncertainty in the energy scale. With a more 
precise calibration method, which will be available in the real experimental conditions, it 
is reasonable to think that the absolute energy scale will be known to the 1% level, hence 
reducing the systematic error of p to 0.04. Conversely, with the calibration method used 
in this analysis, the error of the absolute energy scale is dominated by the systematics 
on the determination of the length of the muon track segments. Therefore, no further 
improvement of the result can be foreseen by using a larger data sample. The final result 
is: 

p = 0.72 ± 0.06 (stat.) ± 0.08 (syst.) (4) 
which is compatible with the V-A value. 



5.2 Energy resolution 



The high level of agreement between the experimental and the generated muon decay 
spectra allows us to estimate the detector's energy resolution using the simulated sample. 
We have evaluated separately the contributions to the final energy resolution coming from 
the following sources: 

• Electronic noise: This contribution depends on the S/N ratio. For non-correlated noise 
and constant S/N ratio (which is a good approximation of the real experimental con- 
ditions) this contribution is expected to vanish at high energies (i.e. for high numbers 
of wires) . 

• Reconstruction effects: This contribution is intrinsic to the signal extraction method 
we use, and independent of the energy. The main contribution comes from hits whose 
shape is not well reproduced by the fitting function. Such a case may arise for tracks 
with small angle with respect to the drift direction, or unresolved close hits [15]. Other 
minor contributions are due to undetected hits or fake hits created during the automatic 
reconstruction procedure. 

• Calorimetric calibration: As already mentioned in Section 4.2, an extra contribution 
to the energy resolution is expected when correcting for the electron attachment to 
impurities using an average value of the drift electron lifetime. This contribution is also 
energy independent and has been estimated in Section 4.2. In the case of this particular 
analysis, we have also considered an extra contribution to the energy resolution coming 
from the determination of to for the different events. 

The first two contributions have been evaluated with MC samples generated in the 
appropriate conditions (see Figure 7). In both cases, the energy resolution is defined 
^ (-^meas ~ -^Mc)/-^MC) where E'mc ^-nd -E^eas ^^^i respectively, the MC generated and 
measured energy of the electron track, disregarding bremsstrahlung losses (i.e. only the 
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Fig. 7. Energy resolution as a function of the electron energy. 

fraction of the energy lost by the electron by ionization of Ar atoms is compared). For 
every bin, the central value and error of the resolution are obtained from a Gaussian 
fit, with errors of a statistical nature. 



In order to evaluate the energy resolution expected from the reconstruction effects, 
we have generated the MC sample with no electronic noise, and passed it through the 
whole reconstruction chain. In such a case, we obtain a constant resolution of (1.97 ± 
0.05)% (see Figure 7). Conversely, the joint effect of the electronic noise and the recon- 
struction is evaluated with the MC sample generated including the electronic noise. A fit to 
the points is performed by the function {E'^^^^ — E^j^q) / EIj^q = a/ ^ E [MeV] ® b, where the 
best fit is obtained for (errors are of statistical origin) a = (11 ± 1)% and b = (2.5 ±0.3)% 
(see Figure 7) . This result is compatible with the assumption of a vanishing contribution 
of the electronic noise at high energies. 



6 Conclusions 



We have performed the first physics measurement with the ICARUS LAr TPC 
detection technique: the determination of the Michel p parameter from the detailed study 
of the muon decay spectrum using the stopping muon event sample from the ICARUS 
T600 detector test run with cosmic rays. We obtain p = 0.72 ± 0.06 (stat.) ± 0.08 (syst.), 
in agreement with the SM value p = 0.75. This measurement involves the exploitation of 
both spatial and calorimetric reconstruction capabihties of the detector. Therefore, the 
obtained result constitutes a proof of the maturity of the detection technique to produce 
high quality physics results, in view of the operation of the detector in the Gran Sasso 
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underground laboratory. We have also estimated the energy resolution including detector 
and reconstruction effects from a sample of MC simulated Michel electrons, obtaining 



(-^meas ~-^Mc) / ^lic Compatible with 11%/ ^JE [MeV] ®2%, where the raise at low energies 
is due to the electronic noise, and the constant term arises from reconstruction effects. 
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